Hydrodynamic attraction of swimming microorganisms by surfaces 
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Cells swimming in confined environments are attracted by surfaces. We measure the steady-state 
distribution of smooth-swimming bacteria (Escherichia colt) between two glass plates. In agreement 
with earlier studies, we find a strong increase of the cell concentration at the boundaries. We 
demonstrate theoretically that hydrodynamic interactions of the swimming cells with solid surfaces 
lead to their re-orientation in the direction parallel to the surfaces, as well as their attraction by the 
closest wall. A model is derived for the steady-state distribution of swimming cells, which compares 
favorably with our measurements. We exploit our data to estimate the flagellar propulsive force in 
swimming E. coli. 
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The majority of swimming microorganisms involved in 
human functions and diseases are found in geometrically 
confined environments. Spermatozoa in the female repro- 
ductive tract swim in constricted domains []]]. Bacteria 
make their way through host cells and tissues @, and 
aggregate in antibiotic-resistant biofilms on surfaces 

Despite the ubiquitous nature of biological motility 
near surfaces, not much is known about the physical 
consequences of locomotion in a confined environment 
0, HI- Perhaps the simplest observed effect of locomo- 
tion near walls is the accumulation of swimming cells on 
surfaces. In 1963, Rothschild measured the distribution 
of bull spermatozoa swimming between two glass plates 
(separation, 200 jiim). The cell distribution was nonuni- 
form, with a constant density in the center strongly in- 
creasing near the walls [6j. Similar results were later ob- 
tained for human spermatozoa in glass tubes Q ■ Further 
studies for animal spermatozoa pointed out the possible 
important of three-dimensional effects [H, Numeri- 
cal simulations of model cells with two-dimensional beat 
patterns [l(| supported an explanation in terms of cell- 
surface hydrodynamic interactions, a scenario confirmed 
by recent computations for suspensions of simplified low- 
Reynolds number swimmers More recent work fo- 
cused on the change in swimming kinematics near solid 
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In this paper, we study the attraction of swimming 
bacteria by solid surfaces. We measure the distribution 
of non-tumbling E. coli fl6| cells swimming between two 
glass plates in a density-matched fluid, and obtain re- 
sults qualitatively similar to that of Rothschild We 
demonstrate theoretically that the origin for the cell pro- 
file is purely hydrodynamical. Using physical arguments 
based on long-range hydrodynamics interactions between 
swimming cells and surfaces, we show that these interac- 



tions induce a reorientation of the cells in the direction 
parallel to the surface, independently of their initial con- 
dition (position, orientation), and the subsequent attrac- 
tion of the cells by the closest wall. Our model allows 
us to predict the resulting steady-state cell distribution, 
and is exploited to obtain an estimate for the flagellar 
propulsive force in swimming E. coli. 

Our experimental procedure is illustrated in Fig. [1] E. 
coli (smooth-swimming strain HCB-437 (l7| ) is grown to 
mid-exponential phase in T-broth (1% Tryptone, 0.5% 
NaCl), washed three times by centrifugation (2200g for 
8 min), and then resuspended in a motility medium (10 
mM potassium phosphate, pH 7.0, 0.1 mM EDTA). PVP- 
40 (polyvinylpyrrolidone) is added (0.005%) to prevent 
adsorption of cells to glass, and the final suspension is 
combined with Percoll (2:3 ratio) to match the medium 
and cell buoyant densities [2^| . A droplet of the cell mix- 
ture is deposited between two glass coverslips, previously 
cleaned in a mixture of ethanol saturated with potassium 
hydroxide, rinsed with ultra-pure filtered water, and al- 
lowed to air dry. The coverslips are separated by a dis- 
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FIG. 1: Representation of the experimental procedure, (a): 
Smooth-swimming E. coli cells are mixed with a density- 
matched fluid; (b): The cell mixture is deposited between two 
glass plates (separation distance H); (c): The distribution of 
swimming cells is imaged as a function of the distance, y, from 
the lower surface; (d): Example of image obtained from data 
acquisition in the first layer above the glass surface. 
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tance H, controlled by layers of other coverslips (#1.5) 
and verified by caliper measurement. A phase-contrast 
microscope (Nikon Optiphot-2) using 600x magnification 
(depth of field, 4.3 /mi) and equipped with a shuttered 
CCD video camera (Marshall Electronics V1070) set for 
an exposure of 1 ms/framc is used to image the popu- 
lation of swimming cells. The video signal is sent to a 
MacG4 equipped with an LG-3 frame grabber (Scion Im- 
age Corp.) and ImageJ software (NIH, Bethesda, MD). 
We capture 2-second movies at 20 frames per second and 
measure the number of swimming cells by counting cells 
swimming at speed larger than 1 body length per sec- 
ond. We start 5 /im above the lower glass surface; we 
then bring the plane of focus up 10 /im, and repeat the 
measurement until we reach within 5 //m of the upper 
glass surface. Experiments are then repeated with other 
cell samples and sets of coverslips. 

In our protocol, two parameters can be varied: the 
distance, H, between the two coverslips (we chose H — 
100 /im or 200 /im) and the cell density of the final mix- 
ture, i.e., the size of the overall cell population (when 
H = 100 /im, we performed additional experiments dou- 
bling the number of cells). The experimental results are 
shown in Fig. [21 vertical errors bars represent statistics 
on ten different experiments, and horizontal error bars 
the depth of field. As in Ref. 0, we find that the cell 
profile peaks strongly near the walls, with a nearly con- 
stant cell density about 20 /im away from the walls; this 
is the main experimental result of this paper. 

We now turn to the physical understanding of the at- 
traction phenomenon. In order to provide a complete 
physical picture, we need to identify the mechanism re- 
sponsible for the non-uniform cell distribution, and pre- 
dict the steady-state profile observed experimentally. 

The physical mechanism for the attraction is the hy- 
drodynamic interactions between swimming cells and 
surfaces [HI El • The flow around most flagellated swim- 
ming organisms, including spermatozoa cells or bacteria 
such as E. coli, is well approximated by a force-dipole 
(stresslet) [l8j |: the flagellar motion provides the propul- 
sive force which is opposed by the drag on both the 
cell body and flagella, corresponding to a force-dipole in 
which both the flagella and the body act on the fluid in 
the direction away from the cell (represented in Fig. [3^ 
by two arrows pointing in opposite directions) . The fluid 

velocity is given by u = (-1 + 3^-) r, where 

p > is the dipole strength, e the swimming direction, 
r\ the viscosity, and r the distance to the dipole; this far- 
field model is valid for distances larger than the length, 
L, of the swimming cells (body plus flagella), an approx- 
imation that we will make in this paper. 

Near a wall, the flow field induced by the cell is a su- 
perposition of that due to the force-dipole, plus any im- 
age flow field, located on the other side of the surface, 
and necessary to enforce the correct surface boundary 
condition (similar to the method of images in electro- 
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FIG. 2: Experimental data: number of swimming cells, n, as 
a function of the distance to the bottom coverslip, y, when 
the distance between the surfaces is H — 100 /im (top) and 
H — 200 /im (bottom). The lines are fit to the data with 
the model of Eq. ([6]) with no = 1.5 and L± = 34.8 (im (top, 
solid line), no = 0.3 and L± — 59.1 um (top, dashed line) and 
no = 3.9 and L± = 26.4 /im (bottom, solid line). 



statics, only here the image is a vector field). The image 
system for a force-dipole parallel to a no-slip surface is 
known [l|| (force-dipole, force-quadrupole and source- 
quadrupole) and is found to induce, at the location of 
the dipole, a velocity component towards the solid sur- 
face of order ~ p/rjU 2 where y the distance to the sur- 
face (Fig. [3)3); this wall-induced flow is the reason for 
the attraction fllj| . To gain physical intuition, it is eas- 
ier to picture a dipole near a free surface; in that case 
the image system is an equal dipole on the other side of 
the surface, and two parallel dipoles attract each other. 
Physically, this wall-induced velocity is reminiscent of 
the shear-induced migration of deformable bodies and 
polymers away from boundaries [13, HD, 22, Hj3], with 



the difference that in these cases, the dipoles arise from 
shear-induced deformation, and therefore they have the 
opposite signs (p < 0). 

Although the aforementioned dipole-dipole attraction 
contains the essential physical picture, some further con- 
siderations are required. When the force-dipole is not 
aligned with the nearby surface, the attraction can be- 
come a repulsion. Indeed, if we denote by 9 the ori- 
entation of the dipole w.r.t. the vertical direction, the 
induced velocity in the direction away from the no-slip 
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FIG. 3: Attraction of microorganisms to solid surfaces, (a) 
The flow field around a swimming cell is well approximated 
by a force-dipole of strength p > 0, represented by two ar- 
rows pointing in opposite directions; (b) Hydrodynamic at- 
traction of a force-dipole by a no-slip surface due to the image 
system on the other side of the surface (force-dipole, force- 
quadrupole, source-quadrupole) ; (c) Notations for the model. 



wall is given by (Fig. [3b) 
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Notably, if the angle is small enough (6 < arccos 1 / v3) , 
the induced velocity changes sign and turns into a re- 



pulsion (lly > 0) 



Moreover, the average of the 



wall-induced attraction over a whole population of ran- 
domly oriented microorganisms is exactly equal to zero, 



Since, for E. coli, p > and 7 > 1 (7 ~ 2 for the cell 
body, and 7 3> 1 for the flagellar bundle), Eq. ([3]) shows 
that f2 z has always the same sign as — cos 8 sin 9. When 
< 6 < 7r/2, the rotation rate is negative, and brings the 
cell in the direction parallel to the surface; when tt/2 < 
9 < 7r, the rotation rate is positive, and also aligns the 
swimming cell parallel to the surface. Consequently all 
swimmers are reoriented in the direction parallel to the 
surface by hydrodynamic interactions [3lj |. 

On what time scale do the cells reorient? Cells not 
initially parallel to the surfaces first swim towards one 
surface, on a time scale t\j ~ H/U, where U is their 
swimming speed. When they reach a distance y ~ L 
from the surface, reorientation takes place on a time 
scale tq ~ Sl^ 1 ~ r]L 3 /p. From a scaling standpoint, 
the dipole strength, p, is on the order of the drag (or 
thrust) force on the organism times the typical cell size, 
p ~ rjUL 2 , and therefore we have tq ~ L/U. For E. coli 
cells of size L ss 5 — 10 /im swimming at U ~ 20 /im s _1 
[iH, the reorientation occurs in a matter of seconds. 
By comparison, the time scale tr for reorientation by ro- 
tational Brownian motion is tr ~ , where Dr is the 
rotational diffusivity, Dr ~ ksT/riL 3 [kg is Boltzmann's 
constant and T the temperature). For E. coli cells, we 



J u y su\9&9 = 0. Consequently, any asymmetry in the fi- expect t r ~ 10 2 s, and therefore the swimming direc- 



nal distribution of swimming cells due to this mechanism 
alone will reflect the non-uniformity of the initial distri- 
bution of cell orientations, and is therefore expected to 
be small for large populations. 

Here, we propose that hydrodynamic interactions with 
the surfaces provide the physical mechanism to modify 
the orientations of the cells. Let us consider the complete 
flow field due to a swimming cell near a solid surface. Not 
only does the image system for the force-dipole induces 
a local attractive/repulsive velocity for the cells, it also 
contains (in general) non-zero velocity gradients; these 
gradients are responsible for the rotation of the cells. 
Modeling the swimming cells as force-free and torque- 
free prolate spheroids of aspect ratio 7, their rotation 
rate, fi, is given by 



tions of the cells are dominated by the (deterministic) 
hydrodynamic mechanism outlined above. 

With cells swimming parallel to the surfaces, the 
steady-state cell probability distribution across the cham- 
ber, n(y), is described by a balance between the advec- 
tion by the velocity field in Eq.([T]) and cell diffusion (with 
diffusivity, D±), so the cell conservation equation is writ- 
ten as 
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(4) 



For cells parallel to two surfaces separated by a distance 
i?, the attractive velocity is written as 
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where u> = V x u and E = ^(Vu + Vu T ) denote, re- 
spectively, the vorticity and the rate-of-strain of the flow 
field due to the image system 24] . Eq. @ states that the 
cells are rotating at one half the local value of the vor- 
ticity generated by the image system, plus an additional 
term that depends on the aspect ratio of the cell and the 
straining flow component of the image velocity field. 

Using the notations defined in Fig. [3l we evaluate the 
component of the rotation rate in the direction 2, parallel 
to the surface and perpendicular to the dipole, and obtain 



plus terms of higher order in 1/y and which are 

due to reflections of each ima ge s ystem on the opposite 
surfaces and are neglected [3,123] ■ Integration of Eq. (j4|) 
using Eq. ([5]) and assuming that dn/dy = halfway be- 
tween the two plates leads to an analytical model for the 
steady-state concentration profile as 
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The comparison between our data and this model are 
presented in Fig. [21 We fit the two parameters in Eq. © , 
no and L, to the data points located further away from 
the closest wall than h ~ L w 10 urn, below which our 
point-dipole model is no longer valid [321 ] . The agreement 
between the theory and our experimental results is good. 
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Close to the wall, we have (y, H — y) <C L±, and the 
model over-predicts the cell density. This could be reg- 
ularized by modeling near- wall hydrodynamics [H, |28| . 
and including other cell-wall interactions (intermolecular 
and screened electrostatics), but would not modify our 
far-field results. For the data with H = 100 fim, the best 
fits with our model (in a least-square sense) are obtained 
for L = 48.3 ± 15 /im; when H — 200 /im, agreement is 
obtained for L = 26.9 ± 4 /um. 

We now exploit our results to estimate the dipole, p, 
for smooth-swimming E. coli cells. From Eq. ([6]), we see 
that for a fluid of known viscosity, ?y, one only needs the 
value D± and the estimate of L, to obtain the estimate 
p = 6Attt]D±L/3. Fitting the measurements to our model 
leads to L = 20 — 76 /im. Since the cells swim parallel to 
the surfaces, the diffusion coefficient, D±, is the Brown- 
ian diffusivity for bacteria in the direction perpendicular 
to their swimming direction [33|; in water [3J] at 20° C, 
we estimate D± rs 0.1 — 0.2 /j,m 2 /s, leading to the es- 
timate p w 0.1 — 1 pN /im. The distance, £, between 
drag- and thrust-producing units on the swimming cell 
is on the order of the size of the cell body, which dis- 
places the center of drag ahead of the center of thrust, 
and £ « 1 /xm. We therefore estimate from our measure- 
ments a thrust force, / = p/l w 0.1 — 1 pN. This agrees 
with recent experiments for swimming E. coli reporting 
flagella thrust forces / w 0.57 pN [if and / « 0.41 pN 

In conclusion, we have studied experimentally and the- 
oretically the attraction of swimming micro-organisms by 
surfaces. Experimental data with smooth-swimming bac- 
teria display a strong attraction by solid surfaces, which 
we have rationalized as follows: hydrodynamic interac- 
tions with surfaces result in a reorientation of the swim- 
ming cells in the direction parallel to the surfaces, and 
an attraction of the aligned cells by the nearest wall. 
We have also shown how to exploit the measurement of 
steady-state population profile to estimate the flagellar 
propulsive force of the swimming cells, a simple method 
to estimate motile propulsive forces. As an extension, we 
note that some swimming cells, such as the algae Chlamy- 
domonas, are not pushed from the back by their flagella 
as in E. coli but are instead pulled from the front. In 
that case, the sign of the dipole is reversed, p < 0: the 
wall-induced rotation rate (Eq. [3|) changes sign, and hy- 
drodynamic interactions reorient the cells in the direc- 
tion perpendicular to the surface. For these cells, hydro- 
dynamic attraction is therefore also expected, but for a 
different reason: they simply crash into the walls. 
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